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Ph' abstract 

00 ■ We review recent developments in the theoretical description of inclusive single- 

ly-} | hadron production at next-to-leading order in the parton model of quantum 

chromodynamics. Fragmentation functions are extracted from fits to data of in- 
clusive pion and kaon production in e + e~ annihilation at different centre-of-mass 
energies. Exploiting sensitivity to the scaling violation, one can simultaneously 
fit the asymptotic scale parameter A so as to obtain an independent determi- 
i nation of the strong coupling constant a s . Owing to the factorization theorem, 

Q\ \ the fragmentation functions only depend on the species of the produced parti- 

cles, but not on the process by which they are produced. This allows one to 
| make absolute theoretical predictions for inclusive pion and kaon production in 

other types of experiments such as photon-photon, photon-hadron, or hadron- 
hadron scattering. Recent data of photoproduction taken by the HI and ZEUS 
Collaborations at DESY HERA nicely agree with such next-to-leading-order 



On 

predictions. 



1. Introduction 



X 1 The La grangian of quantum chromodynamics (QCD) contains quarks and gluons 

as elementary fields. Allowing for these particles to appear as asymptotic states, we 
can evaluate scattering amplitudes perturbatively, in principle with arbitrary preci- 
sion. Of course, this picture needs to be complemented by the principle of confinement 
of colour; experiments detect hadrons rather than quarks and gluons. Nevertheless, 
this simplified computational procedure is very successful in describing the production 
of jets of hadrons at high centre-of-mass (cm.) energies (v^) i n e + e~ annihilation 
and at high transverse momenta (px) in scattering processes. Due to parton-hadron 
duality, clustering partons in the final state according to certain jet definitions yields a 
useful approximation, although this does not account for any details of hadronization. 

On the other hand, experiments are providing us with copious information on 
the inclusive production of single hadrons, which cannot be interpreted along these 
lines. In this case, we need a detailed concept for describing how partons turn into 
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hadrons. In the framework of the QCD-improved parton model, this is achieved 
by introducing fragmentation functions (FF's), D^(x,fi 2 ). The value of D^(x,fj, 2 ) 
corresponds to the probability for the parton a produced at short distance to 
form a jet that includes the hadron h carrying the fraction x of the longitudinal 
momentum of a. Unfortunately, it is not yet understood how the FF's can be derived 
from first principles, in particular for hadrons with masses smaller than or comparable 
to the asymptotic scale parameter, A. However, given their x dependence at some 
scale [A, the evolution with \i may be computed perturbatively in QCD using the 
Altarelli-Parisi0 (AP) equations. This allows us to test QCD quantitatively within 
one experiment observing single hadrons at different values of yfs (in the case of e + e~ 
annihilation) of pt (in the case of scattering). Moreover, the factorization theorem 
guarantees that the functions are independent of the process in which they have 
been determined, and represent a universal property of h. This enables us to make 
quantitative predictions for other types of experiments as well. To summarize, having 
extracted FF's from fits to experimental data, we may test their [i dependence as 
predicted by the AP equations and their universality as postulated by the factorization 
theorem. 

After the pioneering leading-order (LO) analyses of pion, kaon,! and charmed- 
mesoni FF's in the late 70's, there had long been no progress on the theoretical side 
of this field. This may partly be attributed to the advent of general-purpose Monte 
Carlo (MC) event generators based on LO parton-level matrix elements in the early 
80's, which were soon to become very popular in the experimental community. In 
these computer programs, the fragmentation into hadrons is simulated according to 
certain phenomenological model assumptions, e.g., the cluster algorithm in the case 
of HERWIG or the LUND string model in the case of PYTHIA. Although such MC 
packages often lead to satisfactory descriptions of the data, from the theoretical point 
of view, their drawback is that this happens at the expense of introducing a number of 
ad-hoc fine-tuning parameters, which do not originate in the QCD Lagrangian. Fur- 
thermore, in the MC approach, it seems impossible to implement the factorization 
of final-state collinear singularities, which impedes a consistent extension to next- 
to- leading order (NLO). On the contrary, in the QCD-improved parton model, such 
singularities are absorbed into the bare (infinite) FF's so as to render them renormal- 
ized (finite) in a way quite similar to the procedure for the parton density functions 
(PDF's) at the incoming legs. Therefore, a meaningful quantitative test of QCD can 
only be performed in the parton model endowed with FF's at NLO. 

Some time ago, NLO FF sets for 7r°,§ 7r ± , K ± , and rj mesons^ have been con- 
structed through fits to data of e + e~ annihilation generated with HERWIG. By con- 
trast, our procedure, and partly that of Ref. 4, has been to fit to genuine experimental 
e + e~ data. We have introduced LO and NLO FF sets for charged pions and kaonsiS 
and for neutral kaonsl The results presented here are obtained with the up-to-date 
sets of Refs. 7 and 8, which are based on data from SLAC PEP with yfs = 29 GeV 
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and CERN LEP1 with ^J~s = Mz- In contrast to Refs. 4 and 5, Refs. 6-8 also provide 
ready-to-use parameterizations of the /i dependence of the FF's. 

This presentation is organized as follows. In Section 2, we shall report the es- 
sential features of our fitting procedure for charged pions and kaonsB and neutral 
kaons.! In Sections 3 and 4, we shall confront recent data of inclusive pion and kaon 
production collected by the Hl!00 and ZEUS Collaborations^ in jp scattering at 
DESY HERA and by the UA1 Collaboration^ in pp scattering at the CERN SppS 
collider, respectively, with NLO predictions based on our FF's. Our conclusions will 
be summarized in Section 5. 

2. Fit to data of e + e~ annihilation 

The most direct way to obtain information on the FF's of hadrons is to analyze 
their energy spectrum measured in e + e~ annihilation, where the theoretical predic- 
tions are not obscured by additional nonperturbative input, e.g., in the form of PDF's 
for the incoming particles. At NLO in the parton model with nj massless quark 
flavours, the differential cross section of e + e~ — > j/Z — > h + X, normalized to the 
total hadronic cross section cr tot , is given by 
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where x = 2Ef l /y r s is the fraction of the beam energy carried by h, a = g, q±, . . . , q nf , 
jj, is the renormalization scale, and Mh is the factorization scale. The parton-level 
cross sections read 
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where e qi is the effective coupling of qi to the photon and the Z boson including 
propagator adjustments and C a are the NLO corrections ll Here, -P a ^P are the LO 
terms of the timelike a — > b splitting functions, 
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which control the /i evolution of the FF's via the AP equations 
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Analytic expressions for P^J, b may be found, e.g., in Ref. 16. The integro-differential 
equation (|4]) may be solved either via the Mellin transform technique or with brute 
force as it stands. Good agreement is found between these two methods. @ 

It is natural to choose fi = = y/s. This eliminates the terms involving 
in Eq. (0). In e + e~ annihilation, gluon fragmentation occurs only at NLO and be- 
yond. To increase the sensitivity to Dg, one may select longitudinal polarization,^ 
so that the delta function in Eq. (0) does not contribute, or concentrate on gluon- 
tagged three-jet events.00 Then, the respective NLO terms in Eq. (0) constitute the 
Born approximation, and it is desirable to include the next-to-next-to-leading-order 
corrections^ in order to reduce the scale and scheme dependences. Unfortunately, 
these were not yet available when the analyses of Refs. 6-8 were carried out. 

The fitting procedure in Refs. 7 and 8 was as follows. For each of the hadron 
channels h = 7r + + 7r~, K + + K~ , K s + K\ and each parton type a — g,u, d, s, c, b, 
we made the ansatz 

D h a {x,^) = Nx a {l-xf, (5) 
at the respective starting scale 

{V2 GeV if a = g,u, d, s, 

M(r] c ) = 2.98 GeV if a = c, (6) 
M(T) = 9.46 GeV if a = b. 

According to the flavour composition of h, we respectively imposed 

D? (x if j$)=Df (x,iil), 



M) = ^f M). (7) 



The analysis of Ref. 7 was based on charged-pion, charged-kaon, and uniden- 
tified charged-hadron dataMIll Apart from charged pions and kaons, mainly protons 
and antiprotons contribute to the charged-hadron yield. Inspired by Ref. 21, we 
approximated 

da h± . ..da** da K± 

_ = [! + /(,)]_ + _, (8) 

where f(x) = 0.195 — 1.35 (x — 0.35) 2 . In their charged-hadron analysis, the ALEPH 
Collaboration^ distinguished between uds-, c-, and 6-enriched samples. Also exploit- 
ing information on identified gluon jets,00 we were thus able to treat g,u, s,c,b — > 
and g,u, d,c,b — > K ± fragmentation separately. Due to the large gap in y/s between 
PEP@ and LEPl,0'lli'0 we could simultaneously determine A™. Thus, we had a total 

of 2(tt, K) x 5(partons) x 3(N, a, (3) + 1(A^.) = 31 independent fit parameters. These 
turned out to be tightly constrained by our LO and NLO fits. In fact, we obtained 
X 2 /d.o.f. values of 134.4/136 = 0.99 and 125.3/136 = 0.92, respectively. We found 
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Agi = 108 MeV (227 MeV) at LO (NLO), which corresponds to a,(M%) = 0.122 
(0.118). This nicely agrees with the latest LEP value, (0.121 ± 0.003)11 

In our neutral-kaon analysis,! we adopted the values from Ref. 7. Further- 
more, appealing to the flavour blindness of the gluon, we assumed that 

Df M)=J?f M), (9) 

at /i = v^2 GeV. Thus, the number of independent fitparameters was 4(partons) x 
3{N,a,P) = 12. Our combined fit to the MARK iS and ALEPH0 neutral-kaon 
samples yielded xVd.o.f. = 9.9/20 at LO and 8.6/20 at NLO. 



3. Comparison with data of photoproduction in ep scattering 

According to present HERA conditions, E e = 27.5 GeV positrons collide with 
E p = 820 GeV protons in the laboratory frame, so that ^fs = 300 GeV is available in 
the cm. frame. It has become customary to take the rapidity of hadrons travelling 
in the proton direction to be positive. The rapidities measured in the ep laboratory 
and cm. frames are related through 

Vcm. =2/lab- J^TT- ( 10 ) 

In photoproduction, the electron or positron beam acts like a source of quasi-real 
photons, with low virtualities — Q 2 , so that HERA is effectively operated as a 7p 
collider. The appropriate events may be discriminated from deep-inelastic-scattering 
events by electron tagging or anti-tagging. The photon flux is well approximated by 
the Weizsacker- Williams 2 ^ formula, 



F e (x O 2 ) = — 



l + (l-*) 2 ln gk + 2m 2 g fj_ 

m +zm e xi 




(11) 



where x = E 1 /E ei Q 2 min = m 2 e x 2 /(l - x), and Q 2 max = 0.01 GeV 2 (0.02 GeV 2 ) for 
tagged events in the case of HI (ZEUS). The cross section of ep — > h + X emerges 
from the one of ^p — > h + X by convolution with F^. By kinematics, x m j n < x < 1, 
where x m i n = prexp(— y c .m.)/[\/s — pr ex p(2/c.m .)]• In conformity with the HI and 
ZEUS event-selection criteria, we impose 0.3 < x < 0.7 and 0.318 < x < 0.431, 
respectively. 

It is well known that ^p — > h + X proceeds via two distinct mechanisms. The 
photon can interact either directly with the partons originating from the proton (di- 
rect photoproduction) or via its quark and gluon content (resolved photoproduction). 
Both contributions are formally of the same order in the perturbative expansion. 
Leaving aside the proton PDF's, F£, and the FF's, D 1 ^, which represent common 
factors, the LO cross sections are of 0(aa s ) in both cases. In the case of the resolved 



5 



mechanism, this may be understood by observing that the ab — > cd cross sections, 
which are of 0(a 2 ), get dressed by photon PDF's, whose leading terms are of 
the form aln(M 2 /A 2 ) oc a/a s , with M 7 being the corresponding factorization scale. 
Here, a, b, c, d denote quarks and gluons. In fact, the two mechanisms also compete 
with each other numerically. Resolved photoproduction dominates at small pt and 
positive 7/iab, while direct photoproduction wins out at large px and negative y^. 

LO calculations suffer from significant theoretical uncertainties connected with 
the freedom in the choice of the renormalization scale, /i, in a s and the factorization 
scales, M 7 , M p , and M h , in FJ, F£, and D%, respectively. In order to obtain reliable 
predictions, it is indispensable to proceed to NLO. Let us first consider resolved 
photoproduction, which is more involved. Starting out from the well-known LO cross 
section of ■yp — > h + X, one needs to include the NLO corrections, K ab ^ c , to the 
parton-level cross sections, da a b^ c /dt, to substitute the two-loop formula for a s , and 
to endow F£, and D% with NLO evolution. This leads to 
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-K a ^ c (s, t, u, /x 2 , M 2 , M 2 , M 2 ) (l + 



(12) 



where a, b, c = g,q u . . .,q nf , s = (p a + p b ) 2 , t = (p a - p c ) 2 , and u = (p b - p c ) 2 . 
The parton momenta are related to the photon, proton, and hadron momenta by 
p a = x 7 p 7 , p& = XpPp, and p c = Ph/xh- The K a ^ c functions may be found in Ref. 26 
for M 7 = M p . This restriction was relaxed in Ref. 27. 

The NLO cross section of direct photoproduction emerges from Eq. ([12|) by sub- 
stituting F a y (x 7 , M 2 ) = 5(1 — x 7 ), replacing da a b^c/dt and K a ^ c with da^^c/dt and 
fT 7 b_> c , respectively, and omitting the sum over a. The K lb ^ c functions were first 
derived in Ref. 28 setting M 7 = M p = M^ and taking the spin-average for incoming 
photons and gluons to be 1/2. In Ref. 27, the scales were disentangled and the spin- 
average convention was converted to the MS scheme, i.e., to be l/(n — 2), with n 
being the dimensionality of space-time. Analytic expressions for the K lb ^ c functions 
are listed in Ref. 29. 

The theoretical predictions presented here are calculated at NLO in the MS scheme 

(5) 

with rif — 5 quark flavours and = 202 MeV. Unless otherwise stated, we use the 
CTEQ4M0 proton PDF's, the GRV§ photon PDF's, and the BKK FF's for charged 
hadrons@ and neutral kaons.i We set /i = M 7 = M p = Mh = £pr, with £ = 1/2, 1, 2, 
in order to estimate the theoretical uncertainty. 

In Fig. H|, we compare the ZEUS0 and preliminary Hli data on ep — > h^ 1 + X via 
photoproduction with the corresponding NLO predictions. We find good agreement 
as for both normalization and shape. In the upper p? range, our central prediction 
slightly overshoots the ZEUS data (see Fig. |l]a), while it tends to be a tiny bit below 
the centres of the HI data points, but well within their errors bars (see Fig. [l|b). 
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Fig. 1. Comparison of ZEU9I3 and preliminary Hlu data on ep — -> ft. 1 * 1 + X via photoproduction with 
NLO predictions. 



The study of the y^ spectrum in Fig. |I]c nicely illustrates the interplay of direct and 
resolved photoproduction. The direct-photon (resolved-photon) contribution peaks 
at negative (positive) yi a b- At NLO, both contributions strongly depend on the fac- 
torization scheme and scale associated with the incoming photon leg, while their 
sum represents a meaningful physical observable. In the MS scheme with £ = 1, 
the resolved-photon contribution clearly dominates for a minimum-p^ cut as low as 
2 GeV. This offers the opportunity to probe the photon PDF's. In particular, is 
only feebly constrained experimentally. In the forward direction, at yu^ > 1, it makes 
up more than 50% of the cross section (see Fig. |l|d). Thus, a dedicated experimental 
study in the forward direction could help to pin down F£ and to distinguish between 
the various available photon PDF sets. Figure |]d also shows the predictions for the 
acfgpH set with massless charm quark and for the updated GsH set. Unfortu- 
nately, the minimum-pr cut in the data is still too low to allow for a sufficiently 
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precise theoretical description within the Q CD-improved parton model. 

The H100 data on ep — > K°/K° + X via photoproduction are in good agreement 
with the corresponding NLO predictions (see Fig. |2]). In Fig. 0b, the vertical error bars 
include the statistical errors quoted in Table 4.11 of Ref. 11 and an overall systematic 
uncertainty o 10% added in quadrature. Strictly speaking, the minimum-pr cut of 
2 GeV is too low for the parton-model prediction to be usefully precise. Nevertheless, 
Fig. Qb nicely exposes the potential of the yi a b distribution to probe FJ and . In 
view of the general circumstance that the gluon FF's are less tightly constrained by 
the e + e~ data than the quark FF's and the particular assumption underlying Eq. (|9]), 
an independent experimental determination of would be especially desirable. 

4. Comparison with data of pp scattering 

The high-statistics data on pp — > /r* 1 + X and pp — > Kg + X recently published 
by the UA1 Collaboration@0 offer yet another opportunity to test the universality 
of the FF's predicted by the factorization theorem in a nontrivial way (see Fig. |3J). 
At low pt, the bulks of the cross sections are due to the gluon FF's, which makes pp 
scattering complementary to e + e~ annihilation]! The charged-hadron data are spread 
over a wide range in p?, way up to px = 25 GeV, and thus carry intrinsic information 
on the scaling violation of fragmentation. This is nicely illustrated by the dotted line 
in Fig. §a, which emerges from the upper solid line by suspending the AP evolution 
of the FF's and evaluating them at the fixed scale fio = \/2 GeV instead. This leads 
to a significant increase at high p^, by more than a factor of 5 at pt = 25 GeV. The 
experimental data clearly favour the scaling violation encoded in the solid line. 
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5. Conclusions 

The comparative study of inclusive single-hadron production in e + e~, ep, and 
pp collisions allows for a meaningful quantitative test of the Q CD-improved parton 
model and, in particular, of the scaling violation and universality of fragmentation. 
Furthermore, photoproduction experiments at HERA provide useful information on 
the interplay of the direct- and resolved-photon mechanisms. Extensions of these mea- 
surements to higher values of pt and would render an independent determination 
of a s possible and considerably improve our knowledge of F^. 
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